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ABSTRACT 



Aims. We analyze chemical and kinematical properties of about 850 FGK solar neighborhood long-lived dwarfs observed with the 
HARPS high-resolution spectrograph. The stars in the sample have log g >4 dex, 5000 < T^g < 6500 K, and -1.39 < [Fe/H] < 0.55 
dex. The aim of this study is to characterize and explore the kinematics and chemical properties of stellar populations of the Galaxy 
in order to understand their origins and evolution. 

Methods. We apply a purely chemical analysis approach based on the [ff/Fe] vs. [Fe/H] plot to separate Galactic stellar populations 
into the thin disk, thick disk and high-o- metal-rich (ho-mr). Then, we explore the population's stellar orbital eccentricity distributions, 
their correlation with metallicity, and rotational velocity gradients with metallicity in the Galactic disks to provide constraints on the 
various formation models. 

Results. We identified a gap in the [a/Fe] - [Fe/H] plane for the ^-enhanced stars, and by performing a bootstrapped Monte Carlo 
test we obtained a probability higher than 99.99% that this gap is not due to small-number statistics. Our analysis shows a negative 
gradient of the rotational velocity of the thin disk stars with [Fe/H] (-17 km s"' dex"'), and a steep positive gradient for both the thick 
disk and hamr stars with the same magnitude of about -1-42 km s"' dex"' . For the thin disk stars we observed no correlation between 
orbital eccentricities and metallicity, but observed a steep negative gradient for the thick disk and ho-mr stars with practically the same 
magnitude (r; -0.18 dex"'). The correlations observed for the nearby stars (on average 45 pc) using high-precision data in general 
agree well with the results obtained for the SDSS sample of stars located further from the Galactic plane. 

Conclusions. Our results suggest that radial migration played an important role in the formation and evolution of the thin disk. For the 
thick disk stars it is not possible to reach a firm conclusion about their origin. Based on the eccentricity distribution of the thick disk 
stars only their accretion origin can be ruled out, and the heating and migration scenario could explain the positive steep gradient of 
V^ with [Fe/H]. Analyzing the ho-mr stellar population we found that they share properties of both the thin and thick disk population. 
A comparison of the properties of the homr stars with that of the subsample of stars from the N-body/SPH simulation using radial 
migration suggest that they may have originated from the inner Galaxy. Further detailed investigations would help to clarify their 
exact nature and origin. 

Key words, stars: abundances - stars: kinematics and dynamics - Galaxy: disk 



1. Introduction 

The formation and evolution of the Galactic disks is an important 
topic in contemporary astrophysics. The Milky Way (MW) has a 
composite and complex structure with several main components 
(halo, bulge, thin disk, and thick disk) and many stellar streams 
with origin far from being thoroughly understood. 

The main three stellar populations of the MW in the solar 
neighborhood are the thin disk, the thick disk, and the halo, al- 
though most of these stars belong to the thin and thick disks. The 



subdiv ision between the thick disk and thin disk was first identi- 
fied bv lGilmore & Reidl(ll983h . who analyzed the stellar density 
distribution as a function of distance from the Galactic plane. 
These two populations have different kinematics and chemical 
properties. Generally, the thick disk is thought to be composed 
of relatively old stars (e.g. Bensbv et al.ll2003l l2005t iFuhrmannI 



20081: Adibekvan et al.ll201 lb. metal-poor, and cy - enhanced (e.g. 
FuhrmannI 119981 120081: iProchaska et"aD 120001: iFeltzing et all 



2003^. iRedd vetalJ I2006t iHavwoodl l2008t iLee et al.ri2011l: 



* Based on observations collected at the La Silla Paranal Observatory, 
ESO (Chile) with the HARPS spectrograph at the 3.6-m telescope (ESO 
runs ID 72.C-0488, 082.C-0212, and 085.C-0063). 



Adibekyan et alj 1201 2|) that move in Gala ctic orbits with a 
large-scale height a nd long-scale leng t h (e.g. Robin et aD 19961 : 
iBuser et al.l 119991: Ijuric et al.l l2008l iKordopatis et al.l l201lh . 
However, some of the observed trends may depend somewhat on 
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the c riteria applied to d i vide stars into the two disk populations 
(e.g. iFuhrniamil l2008t ISchonrich & BinnevI l2009bl iLee et al.l 
l2011h . Interestingly, this model has been recently challenged. 
Recent analyses of the geometric decompositions of the Galactic 
disk based on the elemental-abundance sel ection of the samp le 
stars yielded strikingly d if ferent result s (see Bow et al1l2012al[bt 
iLiu & van de VenI l2012h . iBovv et alJ (1201 2ah find that mass- 
weighted scale-height distribution smoothly varies and there is 
no thin-thick disk bi-modali ty (i.e., the MW ha s no distinct thick 
disk). We refer the reader to llvezic et al.l(l2012l) and lRix & Bowl 
(1201 3h for more recent reviews of the stellar disk(s) and popula- 
tions. 

While considering that the thick disk has formed only by 
one single mechanism is most probably too semplicistic, sev- 
eral different scenarios have been proposed for the formation 
of thick disk: heating of a pre-existing old thin disk via m inor 
mergers (e.g. lOuinn et al.|[l993l: IVillalobos & Helmill2008l). di- 
rect a ccretion of stars from disrupted satellites (e.g. lAbadi et aP 
2003h. gas accretio n at high redshift and stars formed in situ (e.g. 
Brook etal. 7005), and radial migration of stellar orbits (e.g. 
Schonrich & Binnev 2009a b; Loebman et al. 201 1). 

The proposed scenarios offer very different observation- 
ally testable signatures, hence certain well-defined properties 
of the Galactic disks should be able to distinguish between 
these possible scenarios. In this paper, to provide constraints 
on the suggested various formation models , we explore the 
stellar orbital eccentr icity distributions (e.g. ISales et aP l2009t 
iDi Matteo et al.l 1201 lb . a possi ble correlation b etween the ec- 
centricities and metallicity (see iLee et al.ll201 ih . and rotational 
velo city gradi ents wit h metallicity in the thin and thick disks 
(e.g.lLee et al.i201 h:Kordopatis et al.l201 UlNavarro et al.l20l"lt 



lLiu&vandeVenll2012h . 

For our analyses we used the stellar sample of 1111 long- 
lived FGK dwarf stars from Adibekvan et al. (2012). To separate 
and characterize the different Galactic stellar subsystems, we fo- 
cus on the [a/Fe] ratio (here "ff" refers to the average abundance 
ofMg, Si, andTi). 

The outline of this paper is as follows: In Sect. |2]we intro- 
duce the sample and present the chemical criteria of the disks 
deviation. Results and discussion of comparisons of our results 
with the predictions of different contemporary disk formation 
scenarios follow in Sect. [3] Finally, in Sect.H] we draw our main 
conclusions. 



2. The sample and disk definition 

It is becoming increasingly clear that a dissection of the 
Galactic disks based only on stellar abundances is superior to 
kinematic separation (e.g. Navarro et al. [2011 iLee et al.l 1201 it 
lAdibekvan et al.ll201 UlLiu & van de Venll2012l) . This is because 
chemistry is a relatively more stable property of a star than the 
spatial positions and kinematics. In this analysis, to separate the 
thin- and thick disk stellar components, we used the position of 
the stars in the [c/Fe] - [Fe/H] plane. 



2. 1 . Sample selection and stellar parameters 

Our initial sample comprises 1111 FGK dwarfs with high- 
resolution spectra observed with the HARPS spectrograph 
(iMavor et al.ll2003h at the ESO 3.6-m telescope (La Silla, Chile). 
It is a combination of three HA RPS subsamples: HA RPS- 
1 (IMavor et al.l I200I. HARPS-2 ( iLo Curto et al.l I20T0I) . and 
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Fig. 1. [a/Fe] versus [Fe/H] for the whole sample. Magenta as- 
terisks represent the stars belonging to the halo by their kinemat- 
ics. The black dotted line is the fiducial for division into thin- 
and thick disk populations. The black filled squares refer to the 
thick disk stars, blue triangles to the hamr, and the red dots to 
the thin disk stars. On the top-right corner the high-a and low-a 
separation histogram, after subtracting the separation curve, is 
plotted. 



stars, 376 of which hav e been selected from the CORALIE 
volume-limited samplq^ dUdrv et al.ll2000l) for being non-active 
and low-rotating and the rest are southern confirmed exoplanet 
host stars. HARPS-2 sample was compiled with stars in the so- 
lar ne ighborhood out to 57.5 pc from the Sun (iLo Curto et al.l 
120101) to complete a planet-search survey previously started with 
the CORALIE spectrograph. Then active, high-rotating, known 
binaries and va riable stars \yere ex cluded which lead the sam- 
ple to 582 stars (ISousa et al.ll201 la) . It is worth to note, that for 
the CORALIE sample a color cut-off (B-V < 1.2) was made. 
HARPS-4 sample con sists of 97 star s and was compiled from 
the catalog of Nordstro m et al.l (l2004l) . selecting stars with b - y 
> 0.33 and with ph otometric [Fe/H] between -1.5 and -0.5 dex 
( ISousa et al1l201 la) . We note that there are four stars in com- 
mon between HARPS -1 and HARPS-4, and 14 stars in common 
between the HARPS-2 and HARPS-4 subsamples. 

Precise stellar parameters (Teff, log g, [Fe/H], and ^t) and 
elemental abundances for 12 elements (Na, Mg, Al, Si, Ca, 
Ti, Cr, Ni, Co, Sc, Mn, and V) for all the stars were deter- 
mined in an homogeneous manner Briefly summarizing, the 
equivalent widths of the lines were automatically measured 
with the ARESQ code (Aut omatic Rout i ne for line Equivalent 
widths in stellar Spectra - ISousa et all (|2007|) ). Then the at- 
mospheric parameters and elemental abundances were deter- 
mined using a local thermodynamic equilibrium (LTE) anal- 
ysis relative to the Sun with the 2010 revised version of the 



HARPS-4 dSantos et alJ[20Tlh . HARPS- 1 is composed of 451 



' The Umit distance for the stars in this sample depends on the stellar 
spectral type, i.e., a color-dependent distance cutoff was applied. 

^ The ARES code can be downloaded at 

|http://www.astro.up.pt/sousasag/ares| 
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Table 1. The average values of the Ulsr, ^lsr, and Wlsr velocity components and their standard deviations for the stellar groups, 
along with the number of stars and their percentage. The values obtained from the Gaussian fiting of the data presneted in brackets. 



Stellar groups Ulsr ctu 



O-y 



W, 



o"w 



[Fe/H] 



N 



Percentage 



Thick 


-5(-5) 


66(56) 


-52(-53) 


39(42) 


-8(-22) 


44(36) 


-0.60±0.23(-0.58) 


84 


9.94±1.03 


hamr 


-7(-3) 


42(39) 


-23(-18) 


26(28) 


0(-2) 


21(20) 


0.03±0. 12(0.03) 


58 


6.85±0.86 


Thick+hffmr 


-6(-5) 


58(45) 


-40(-35) 


37(41) 


-4(-10) 


36(32) 


-0.34+0.37 


142 


16.78+1.28 


Thin 


0(-3) 


37(37) 


-9(-6) 


22(23) 


0(-l) 


18(16) 


-0.06±0.22(-0.05) 


692 


81.79+1.32 


Halo 


-39 


154 


-198 


60 


-27 


73 


-0.82±0.23 


12 


1.42+0.40 



spectral synthesis code MOOGu (ISnedenI Il973h and a grid 
of Ku rucz ATLAS9 plane-parallel model atmospheres (iKurucJ 
Il993h . The refer ence abundances used in th e abundance analysis 
w ere taken from Anders & Greyesse (Il989| . We re fer the reader 
to lAdibekvan et al.1 (l2012h and lSousa et al.1 (l2008h for more de- 
tails. 

In order to assemble a sample of stars with the most reli- 
able and precise abundance determinations, we decided to estab- 
lish a cutoff in temperature: we excluded stars with Teff < 5000 
K and T^ff > 6500 K. In those temperature regions the errors, 
both in stellar parameters and rX/H] ab u ndances, are higher (see 
Neves et al.ll2009l: ISousa et al]l20I lallbl: lAdibekvan etani20I2t 



Tsantaki et alj|2013l for more details). At the end of the estab- 



lishing cutoff on temperature we finished with 869 stars. Since 
we are interested in long-lived dwarfs we excluded stars with 
log g <4- dex (23 stars), leading to a final sample of 846 stars. 

Our final sample spans the metallicity range -1.39 < [Fe/H] 
< 0.55 dex, although there are only nine stars with [Fe/H] < -0.9 
dex, and two stars with [Fe/H] > 0.4 dex. The typical relative 
unc ertainties in the metallic ity and [a/Fe] are of about 0.03 dex 
(see lAdibekvan et alj|2012l) . 



2.2. Dissecting ttie Galactic disks 

As mentioned above, the thin- and thick disk stars are different in 
their a content at a given metallicity ([Fe/H]). Here we use this 
dichotomy in the chemical evolution to separate different stellar 
po pulations. 

lAdibekvan et al.l (l20Ilh showed that the stars of this sample 
fall into two populations, clearly separated in terms of [a/Fe] up 
to super-solar metallicities. In turn, high-ff stars were also sepa- 
rated into two families with a gap in both [a/Fe] and metallicity 
([a/Fe] * 0. 17 dex and [Fe/ H] » -0.2 dex) dist r ibutio ns. Here we 
use the same technique as in lAdibekvan et al] ( |20II|) to separate 
the stellar group^. 

The [a/Fe] versus [Fe/H] separation plot for the sample stars 
is depicted in Fig. [Tj The histogram (right-top comer) shows 
that the separation between the stellar families (low-a and high- 
a) is clear. Many recent papers separate the thin and thick disk 
stars by the [a/Fe] "ga p" at a 0.2 - 0.3 dex (see iNavarro et al.l 
I2OIII: iLee et al] l20Ilh and other papers also report the exis- 
tence of the "knee" in [a/Fe] trends for the thick disk stars 
(kinematic ally selected) when [Fe/H] reaches to ^ -0.3 dex (e.g. 
iFeltzing et al. 2003; Bensbv et al. 2007). A similar gap in [a/Fe] 
at X 0.2 dex and similar downturn in [Fe/H] at ^ -0.3 dex are 
seen in Fig.[Tl and as already noted there is also a gap between 
metal-poor and metal-rich a-enhanced stars (the gap in [Fe/H] 
can be seen in top-left panel of Fig. |2]i. We note that there is no 



^ The source code of MOOG can be downloaded at 
|http: //www.as. utexas.edu/~chris/moog. html 

We note that in the current analysis we established a more conser- 
vative temperature cutoff than in our previous work. 



[Fe/H] gap found for the thin disk stars. The presence of a metal- 
rich, old thick disk is already no ticed in the spect r oscop ic stud- 
ies of iFeltzing & Bensbvl ( 120081) and lCasagrandS ( |20I2|) . and a 
m etal-rich g-rich popul ation has also been recently confirmed 
bv lGazzanoetaI.N2()I3h . 

In the [Fe/H] region from -0.28 to -0.18 dex there are only 
three stars with enhanced a abundances, while at lower and 
higher metallicities the number of a-enhanced stars is higher (9 
and 13, respectively). In order to analyze whether this lack of 
high-a stars in the mentioned region is statistically significant or 
is due to a small-number statistics, we performed a simple boot- 
strapped Monte Carlo test. We randomly drew 119 stars, which 
is the number of high-a stars in the metallicity region -0.7 < 
[Fe/H] < 0.3 dex, from our sample in the same metallicity region 
(independent on their a-content). Then we counted the number 
of stars in the interested metallicity region: from -0.28 to -0.18 
dex. We also varied this interval by 0.05 dex ([-0.33; -0.23] to 
[-0.23; -0.13]). The entire process was repeated 10^ times. This 
test indicates that such a small amount of stars in the -0.28 to - 
0. 18 dex metallicity region (or in the above mentioned intervals) 
can be obtained by chance with less than 0.01% probability for 
four (or less) stars, and 0.05% probability for five (or less) stars. 
In our simulation we had only one case where the mentioned 
metallicity interval contained three or less stars (0.001%). This 
Monte Carlo test indicates that it is very unlikely that the ob- 
served gap between metal-poor and metal-rich a-enhanced stars 
is due to small-number statistics and therefore the separation is 
considered real. 

As mentioned before, our sample consists of three HARPS 
subsamples, one of which (HARPS-4) was built to include only 
metal-poor stars. Having this in mind, one might suspect that the 
observed gap could be a result of over-sampling of a thick disk 
metallicity distribution. However this can not be the case, be- 
cause our thick disk sample includes only 32 stars from HARPS - 
4, and almost all of them have metallicities below -0.5 dex (there 
is only one star with [Fe/H] > -0.4 dex, and two stars with -0.5 
< [Fe/H] < -0.4 dex). 

The magenta asterisks in Fig. [Tj re fer to stars belonging to 
halo selected by their kinematics (see lAdibekvan et al.ll201 ll) . 
As is clearly evident, halo stars are divided into two high-a and 
lo w-a groups. This dichotom y confirms the results already found 
bv lNissen & Schusteij(l2010l) . 

2.3. Cliaracterisation of tlie Galactic components 

Guided by the chemical separation, we can characterize the stel- 
lar families in terms of their kinematics and metallicity. Fig. |2] 
shows the distributions of the C/lsr, Vlsr, and Wlsr space ve- 
locity components relative to the local standard of rest (LSR) and 
iron abundance for our stellar groups (for details of the compu- 
tation of the space velo city components we refer the reader to 
lAdibekvan et al.l (|2012|) ). In the plot we also separately present 
the distributions of the thick disk + hamr sample stars because 
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[Fe/H] 





-200 -150 




Fig. 2. Distributions of the Galactic space velocities and metallicity of the stars from the thick disk (black dotted), hamr (blue 
dashed), the thin disk (red solid), and the thick disk -i- hornr (green dotted-dashed). Gaussian fits of the data are also presented. The 
mean of gaussian fits is also pictured by vertical lines on the top of each distribution. Note that the distributions of different groups 
were set lower/higher for the sake of clarity. 

Table 2. Observed gradients, correlation coefficients of V^ and eccentricity with metallicity, the number of stars, and significance 
of the correlation for different stellar populations. 



Stellar groups 


dV^/dlFe/H] 
(kms-' dex-') 


■"flV^/fllFe/H] 


N(V^) 


no-(V^) 


de/d[Fe/H] 
(dex-') 


'^fle/fl[Fc/H] 


N(e) 


ncr(e) 


Thick 


41.9±18.1 


0.247 


84 


2.24a- 


-0.184±0.078 


-0.266 


73 


2.24cr 


hamr 


41.4±28.2 


0.190 


58 


1.430- 


-0.185±0.138 


-0.212 


40 


1.30cr 


Thick+ho-mr 


43.9±7.6 


0.435 


142 


5.19cr 


-0.208±0.036 


-0.475 


113 


5.270- 


Thin 


-16.8±3.7 


-0.164 


692 


4.43cr 


-0.023±0.015 


-0.212 


515 


4.81cr 



the link between these two stellar groups is still under debate 
(see Sect. [3]). The average values of the [Fe/H] and L'^lsr, Vlsr, 
and Wlsr velocity components and their standard deviations, 
along with their fractions are presented in Table 1 . We also per- 
formed Gaussian fits of the distributions of C/lsr, ^lsr, VVlsr 
and [Fe/H] for each Galactic component, and present the mean 
values and dispersion in the same table (the values in brackets). 
We did not perform Gaussian fitting for the halo stars because of 
the small number of stars in this group. We also did not fit [Fe/H] 
distribution of thick disk -i- hamr family, because the distribution 
is obviously bimodal and cannot be fitted with one Gaussian. 

From the Fig. [2] and the corresponding table it is obvious 
that the velocity distributions of the hcmr stars are intermedi- 



ate (even more similar to those of the thin dis k stars) between 
thin and thick disk population stars (see also lAdibekvan et al.l 
i2011i) . but with their average velocities and standard deviations 
having intermediate values between the thin- and thick disk fam- 
ilies. The Kolmogorov-Smimov (K-S) statistics predicts 0.52, 
0.01, and 0.02 probabilities (Pks) that hamr and thin disk stars 
come from the same underlying distribution for Ulsr, ^lsr, and 
W'lsr- The same statistical test gives much lower probabiHties 
when we compare the velocity distributions of the ho-mr and 
thick disk stars - Ulsr (Pks ~ 0.38), Vlsr (Pks ~ 3x10-^), and 
W'lsr (Pks ~ 0.005). From the plot it can be noted that the veloc- 
ity distributions of the thick disk stars are not perfectly Gaussian. 
This can be seen also from the Table 1 where the arithmetic av- 
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Fig. 3. The Toomre diagram for the entire sample. The symbols 
are the same as in Fig. 1. 



erages (and nns) and means of the Gaussian fits (and standard 
deviations) are different. The kinematic parameters obtained for 
the thin and thick disks, in general, agree well with those ob- 
tained for kinematical ly selected thin and thick disks in the so- 
lar neighborhood (e.g. iBensbv et ani2003t ISoubiran et"aDl2003t 
[Robin et al. 2003) , but lo cal normalization for the thick disk is 
a bit higher (but see, al so Squbiran et al.ll20()3t iMishenina et al.l 



l2004HKordopatis et al.ll201 Ih . 

For illustrative purposes as a complement to Fig. |2] the dis- 
tribution of stars of our sample in the Toomre diagram is shown 
in Fig. [3] 

Although in Table 1 we present the average values of the ve- 
locities and metallicity obtained for halo stars, we note that these 
values should be considered with caution because the number of 
stars in the sample is very low. 

3. Results and discussion 

As mentioned in the Introduction, several scenarios have been 
proposed for the formation of the Galactic disks which should 
have left different dynamical and chemical imprints. In this sec- 
tion we test the predictions of the mentioned scenarios using our 
sample of local thin- and thick disk FGK dwarfs as identified 
above. At the end of the section we discuss main properties of 
hamr stellar family. 



3. 1 . Rotational velocity and metallicity 

The large dispersion in the metallicities of the thin disk 
stars in the sola r neighborhood, according to the migration- 
based scenario s (Sellwood & Binnev 2002; Roskar et al.l l2008t 
ISchonrich & Binnev. .2009a: .Minchev & Famaev..2010.) . can be 
explained by inward radial migration of metal-poor and out- 
ward migration of metal-rich stars. A consequence of the 
radial migration is a negative c orrelation between rotational 
velocities (Vd>) and metallicity dSchonrich & Binnevl l2009at 
iLoebman et al.l I2OIII) for the thin disk (Iow-[a/Fe]) stars. We 
note that this prediction is not an "exclusive" to radial migra- 
tion, models which assumed inside-out formation can also pro- 
duce this anti-correlation. This trend in the young stars arise 
from epicyclic motions of stars with their birth radii imprinted 
into their metallicities and disappears with time due to com- 







(thin) • avya[Fe/H] = -16.S±3.7 
(hamr) i aV,/a[Fe/H] = 41.4±aa.2 
(thick) ■ aV,/a[Fe/H] . 41.9±1B.1 
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Fig. 4. Galactic rotational velocity (filled symbols) and velocity 
dispersion (open symbols) gradients with metallicity for stars as- 
signed to different stellar populations. The symbols are the same 
as in previous plots. The smaller symbols correspond to the real 
stars and the larger symbols represent the V^ and cry values for 
metallicity bins. Note that the slopes are obtained for the full 
unbinned data. 



plete mixing. Hence, the same models predict small (positive) 
or no correlation between V^ a nd [Fe/H] for stars w ith high- 
[ff/Fe] i.e. thick disk, old stars (ILoebman et al.ll201 Ih . Most of 
the recent observational studies using low-resolution spectra of 
large samples of thin disk stars outside the solar neighbourhood 
(mainly farther than 0.5 kpc from the Galactic disk) show nega- 
tive gradient of V^ with [Fe/H], which ranges from about -20 to 
-40 km s'' dex'' (e.g iLee et al.ll20TTt iLiu & van de Venll20I2l: 
iKordopatis et al.ll201 Ih . The magnitude of the gradie nt also de- 
pends on the distance from the Galactic plane (e.g. Lee et aP 
I2OIII) . Interestingly, a recent study by iNavarro et al.l (l20Ill) 
found little/no correlation between V^ and metallicity for their 
thin disk stars in the solar neighborhood based on a heteroge- 
neous compilation of data from many different sources (their 
sample included also kinematically selected stars). 

To study the degree of radial mixing by exploring the V^- 
[Fe/H] correlation it is important to use a narrow range of galac- 
tocentric radius. Otherwise this correlation will appear naturally 
due to the V^ and metallicity gradients in the Galaxy. Our sample 
includes stars very close to the Sun (on average x 45 pc) and is 
ideal for this search. The top panel of Fig. |4]shows the gradient 
of mean rotation velocitjO with metallicity for different groups 
of stars. The magnitudes of the gradients (the slopes) and the 
correlation coefficients are presented in Table 2. 

The velocity gradient with metallicity obtained for the thin 
disk stars (-17 km s"' dex') is not far f rom the results (-22 km 
s ' dex~') obtained bv lLee et al.l (|20II|) for their thin disk stars 
in the range 0.1 <\ Z |< 0.5 kpc. We remind that this steep gradi- 
ent is in disagreement with the result obtained by Navarro et al] 
(I2OIII) for solar neighborhood thin disk stars. 

Recent studies obtained contradictory results analyzing V^ 
gradient with [Fe/H] for the thick disk s tars. A steep - h40 to 
-1-50 km s ' dex ' gradient obtained by ISpagna et al.l (l20I0h 



and Lee et al.l (l20Ill) u sing the SDSS data was confirmed by 
IKordopatis etal.l(l20 1 lb . who used independent VLT/FLAMES 



' We adopt 8 kpc for the Sun's distance to the Galactic center and 
220 km s"' for the circular velocity of the LSR. 
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Fig, 5. Trends of eccentricities as a function of metallicity for 
stars assigned to different stellar populations. The symbols are 
the same as in previous plots. The smaller symbols correspond 
to the real stars and the larger symbols present the eccentricities 
and standard deviations for each metallicity bin. The slopes are 
again obtained for the unbinned data. 



observations of 700 thick disk s tars. Other studies (e.g. 
Ilvezic et alj [20081: iBond eFall 120101: ILoebman et alj 1201 lb did 
not reveal significant correlation probabl y due to larger e rrors 
in the photometric chemical abundances (iLee et al.ll201lh . We 
note that all the mentioned studies are based on the samples of 
stars located far from the Galactic plane (usually | Z |> 0.5 kpc). 
Fig. |4] and Table 2 indicate a steep and very similar gra- 
dients bot for the thick disk and hamr stars (about +42 km 
s~' dex'). These v alues agree wel l with those obtained in 
ISpagna et all (l2010h : iLee et al.l (l201lh : iKordopatis et"d] (l201lh 
and, as mentioned before, dis agree with the theoretical pre dic- 
tions with radial migr ation by ISchonrich & BinnevI (l2009ab and 
iLoebman et al.l (l201lh . Interes tingly, recent mod el of chemical 
enrichment in the Galaxy by ICurir et aP (1201 2l) assuming an 
"inverse" gradient at high redshifts (z > 3-4) explains the pos- 
itive rotation-metallicity correlation of the old thick disk pop- 
ulation. They showed that using the inside-out formation and 
chemical evolution model of the Galactic disk suggested by 
[Matteucci & Francois (1989|) and [Chiappini et al. (2001), this 
correlation can be established as a result of radial migration and 
heating processes of stars from the inner region of the disk. 



3.2. Stellar orbital eccentricities 

The shape of eccentricity (e) distribution of thick disk stars 
can be used to constrain th e main physical mechanism re- 
spon sible for thier formation (ISaleset al. 2009). Recent studies 
(e.g. lDierickxetal.ll20ld IWilson et all l201lHLee 61^120111 
IKordopatis et alj 1201 Ih Twhich used the ISales et al.l (|2009|) ec- 
centricity test, reached a general conclusion that the eccentricity 
distribution is inconsistent with the accretion scenario, but could 
not distinguish between other published models, i.e. heating, mi- 
gration or merger scenarios. 

To discuss the distribution of the stellar orbital eccentrici- 
ties we cross-m atched our sample with the Geneva-Copenhagen 
Survey sample (ICasagrande et al.ll201 ll) . which provides the or- 
bital parameters of about 650 stars. 
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Fig. 6. Normalized distributions of eccentricities for different 
stellar populations. Designation of the lines kept the same as in 
the previous plots. 



Fig. |5]shows trends of e with [Fe/H] for different stellar fam- 
ilies and the Table 2 presents the observed gradients quantita- 
tively. From Fig.|5]and corresponding table it is obvious that the 
thick disk and hornr stars show very similar trends with the same 
slopes. On the other hand, the trend of the eccentricities for the 
thin disk stars is practically independent on metallicity. Our re- 
sults quantitatively agree well with thos e obtained for th e G-type 
dwarfs from the SDSS/SEGUE survey ( iLee et al.ll20"TTl) . 

Fig.|6|displays the normalized distributions of eccentricities 
for the different subsamples. The e distribution of the thin disk 
stars peak at the first bin (the bin size is 0. 1) and has an average 
value of 0.13. The distribution is also very concentrated around 
small eccentricities, with only 22 stars (4%) having e > 0.3. In 
contrast, the eccentricity distribution of the thick disk stars has 
a wider width and peaks at 0.3 < e < 0.4 bin, with an average 
value of e = 0.31. The eccentricities of the hornr stars have in- 
termediate values between those of the thin- and thick disk stars, 
with an average value a bit closer to that of the thin disk (e = 
0.17). In case of ho-mr stars again very few stars (5 out of 40) 
have e > 0.3 and only one with e > 0.4. Joining the hamr and 
thick disk subsamples we obtain an e distribution with a peak at 
0.2 < e < 0.3 and average value of 0.26. 



A comparison with the model predictions in ISales et al. 



(|2009|) can only rule out the pure accretion scenario (lAbadi et al 
l2003l) . while it does not strongly favour or discard any of the 
other three scenarios. In particular, the shape of the e distribu- 
tion of the thick disk stars does not diff 'er much from the disk 
heating model of lDi Matteo et alJ (1201 ll) . 

Although the hamr stars have smaller e compared to the 
chemically defined thick disk stars, from the Fig. |6]it is diffi- 
cult to conclude that they have different origin or evolution. If 
the hamr stars just represent the metal-rich tail of the thick disk, 
then their smaller eccentricities are a natural consequence of the 
negative gradient of e with [Fe/H] observed for the thick disk 
stars (see Fig.|5]i. 

3.3. Significance of the observed trends 

To evaluate the statistical significance of the inferred correla- 
tions between V^ vs. [Fe/H] and e vs. [Fe/H] we used a bootstrap 
procedure. First, we obtain the zero-centered distribution of the 
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correlation coefficient by randomly bootstrapping (building ran- 
dom samples by shuffling the parameters among the observed 
set of parameters) the observed data pairs 10^ times. We calcu- 
late the correlation coefficient for each of these uncorrelated data 
sets and then the average and standard deviation of these values. 
By assuming a Gaussian distribution for r^y^/^pe/H] i^ae/aiFs/H]^^ 
we can calculate the probability that the Tgy /grpg/jji (^geidWe/nO °f 
the original dataset was obtained by pure chance, but we merely 
present the offset in sigma {ncriV^) (ncr(e)) in Table 2). 

From the table, one can note that the correlations between V^ 
and [Fe/H] for ho-mr and thick disk stars separately are not very 
significant (1.43cr and 2.24cr), even if the slopes are very steep. 
Interestingly, the significance of the correlation increases a lot 
when we consider thick disk+hamr sample (5. 19cr). To check if 
this non-significance in both hamr and thick disk separate sam- 
ples is due to small sizes of the samples or it is due to small 
metallicity ranges, we "merged" the samples in the following 
way. Since the correlation coefficient (and its significance) does 
not change if one add or subtract a constant to the parameters 
(in this case to V^ and [Fe/H]) we shifted the thick disk sample 
(i.e., shifted every point of it) to the locus of hamr stars subtract- 
ing the difference between average values of [Fe/H] and V^ of 
the thick disk and hornr staro- This procedure centers both data 
sets around the same point. Then we calculated the coefficient 
and significance of the correlation using the same bootstrapping 
approach. For the "merged" sample we obtained practically the 
same slope (41.8+14.6) as for the individual samples and corre- 
lation coefficient about 0.234. We found that this correlation is 
significant at about 2.82cr level. If we compare the later value 
with that obtained for the thick disk subsample we can note that 
by increasing the number of objects in the sample by about 70% 
we did not increase significance of the correlation significantly. 
This suggests that the main reason of the obtained "low" signifi- 
cance is probably the small range of [Fe/H] for both independent 
thick disk and especially hamr subsamples. 

We used the same center-shifting approach (as we did for the 
y0-[Fe/H] data pairs) to analyze the significance of the correla- 
tion as a function of the size of the thick disk and hamr sub- 
samples. We found that increasing the size of the subsamples 
(by about 55% for the thick disk stars) does not significantly in- 
crease the significance of the correlations (from 2.24cr to 2.7 Icr 
for the thick disk subsample). This is probably once again due to 
small range of metallicities of the samples. 

In our sample we have few thick disk stars with metallici- 
ties from -1.4 to -0.9 dex and these stars are widely spread in 
this metallicity range. These rather extreme-value stars can in- 
fluence the fitted slope significantly. From Figs. |4]and |5]one 
can note that at lowest metallicities the stars (7 stars with [Fe/H] 
< -0.9 dex) mostly lie above the y0-[Fe/H] correlation line and 
mostly below the e -[Fe/H] co rrelation line. Interestingly, the 
much larger data in iLee et al.1 (12011) suggest similar behavior 
for their most m etal-poor thick disk stars (see Figs. 7 and 9 of 
iLee et al]l201lh . Nevertheless, to check how these data could al- 
ter the slopes we established a cutoff in metallicity at [Fe/H] - 
-1.1, -1, and 0.9 dex. We found that the V'0-[Fe/H] slope changes 
from 48.7+21.9 to 61.5+25.8 km s ' dex"' and the slope of e- 
[Fe/H] varies from -0.26+0.09 to -0.35+0.11 dex"'. This sug- 
gests that decreasing the lower metallicity limit for the thick disk 
the slopes become steeper staying within the one-sigma errors. 



From our tests we can conclude that the observed trends are 
real, but the exact values of the slopes (for the thick disk) should 
be considered with caution. 



3.4. High-a metal-rich stars 

In this subsection we discuss some properties of hamr stellar 
family, which they share with thin- and thick disk population. 

From Fig. |4] it is interesting to see that although the hamr 
and the thin disk stars at [Fe/H] > dex have similar rotation 
velocities, their trends with metallicity are completely different. 
The hamr stars behave like an extension of the thick disk stars 
on the y^-[Fe/H] plane which might mean a similar origin and/or 
evolution. The same similar behavior to the thick disk stars can 
be seen from From Fig. |5] Although the hccmr stars have orbits 
with eccentricities on average more similar to those of the thin 
disk stars (see Fig. |6]l, the e vs. [Fe/H] trend of these stars is the 
same as for the thick disk stars. 

Despite the mentioned similarities to the thick disk, the hamr 
stars also share some properties of the thin disk. From Table 
1 one can see that the dispersion of all the velocity compo- 
nents of hamr stars is very similar to those of the thin disk 
sample. In particular the low dispersion of the vertical veloci- 
ties s uggests a short scale-hei ght, very similar to that of the thin 
disk. lAdibekvan et al.l (1201 ih already showed that the average 
maximum vertical distance (Z^ax) the stars can reach above the 
Galactic plane is about 0.3 kpc for the hamr. For comparison the 
average values of Z^ax are about 0.25 and 1 kpc for the thin disk 
and thick disk, respectively. Adibekvan et al. (201 1) also showed 
that the hamr and thick disk family stars have almost the same 
age, being on av e rage o lder than thin disk stars by about 3 Gyr. 
lAdibekvan et alj ( 1201 2h analyzed the [X/Fe] trends with metal- 
licity of hamr stars for different elements and showed that for 
some elements the trends are different than those of the thick 
disk stars. 

At the same time, hamr stars are more metal-rich than the 
thick disk stars (if they do not just represent the metal-rich tail 
of the thick disk) and they are as metal-rich as the bulge stars 
and have similar enhanced a-element abundances compared to 
Galactic disk stars as recently found for the bulge stars (e.g. 
iFulbright et ai1l2007l iBensbv et al.ll20T3l iNess et al.ll20T3lF l 

The mentioned results hint that hamr stars might have an 
origin in the inner Galaxy (bulge) and have migrated to the sola r 
annulus (see also lAdibekvan et aP 201 1 : Gazzano et al.l l20 131) . 
Very recently, in their simulation lRoskar et al.l(l2012l) identified a 
sub-population with high- [Fe/H] and low- [O/Fe] . Although their 
simulations did not reproduce our data, i.e. separation between 
low- and high-a stars at super-solar metallicities, at first glance 
the higher- [O/Fe] part of that family shares many p roperties of 
the hamr star s. From Fig. 7 of lRoskar et al.l (1201 2l) and Fig. 14 
of Loebman e t al.l (1201 lln we can see that our hamr stars match 
well with the intermediate/old age population of stars with short 
scale-length and scale-height, migrated from the central part of 
the Galaxy (< 2 kpc). 



* We remind thie reader tliat the slopes of the correlation for the two 
samples are very similar. 



' We note that the conclusion that bulge stars are more enhanced in 
a-elements is also based on the Fig. 27 of Bensby et al. 2013, but these 
authors did not conclude that in their paper. 

** We note that in both papers they used the same N-body/SPH simu- 
lation. 
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4. Summary and conclusions 

We selected a sample of about 850 solar neighborhood FGK 
dwarf stars with precise stellar velocity components and chem- 
ical abundances derived from the high-resolution HARPS spec- 
tra. The stars have log g > 4 dex, 5000 < Teff < 6500 K, and -1.4 
< [Fe/H] < 0.55 dex. 

Applying purely chemical approach based on the [ff/Fe] vs. 
[Fe/H] plot, we separate Galactic stellar populations into the thin 
disk, thick disk, and hamr. The later two families are separated 
based on the observed "gaps" in both [Fe/H] and [a/Fe] distribu- 
tions. Performed bootstrapped Monte Carlo test gives a proba- 
bility higher than 99.99% that the observed "gap" is statistically 
significant and is not due to a small-number statistics. 

We characterize the stellar subsamples in terms of their kine- 
matics and metallicity and found in general good agreement with 
earlier determinations of these parameters, although the local 
normalization for the thin disk is about 82 %, which is lower than 
the usually adopted value (but see, als o Soubiran et al.l l2003t 
iMishenina et al.ll2004t iKordopatis et alj|201 Ih . 

Our analysis shows that the rotational velocity of the thin 
disk stars decreases with increasing [Fe/H] (-17 km s"' dex"'), 
while this trend shows an increases with [Fe/H] for both thick 
disk and ho-mr, both with the same magnitude of about +42 km 
s ' dex '. Although the magnitude of the slope for the thick disk 
stars depends on the cutoff established in [Fe/H]. These results, 
obtained for solar neighborhood stars (on average x 45 pc), are 
an independent confirmation of those found for the S PSS sam- 
ple of stars located far from the Galactic plane (e.g. iLee et al.l 
1201 lb . The negative gradient of V^ with [Fe/H] for the thin disk 
qualitatively agrees well with the expectations from the radial 
migration model and the steep positive gradients observed for 
the thick disk andjimnr stars c an be explained by the recent 
simulation bv lCurir et al.l (1201 2l) assuming an inverse chemical 
gradient in the inner early Galaxy. 

For the thin disk stars we observed no correlation between 
orbital eccentricities and metallicity, but observed a steep nega- 
tive gradient for the thick disk and hamr stars with practically 
the same magnitude of -0. 1 85 dex" ' . Again the magnitude of the 
slope for the thick disk stars depends on the cutoff established 
in [Fe/ H]. These results quantitatively agree well with lLee et al.l 
(1201 ll) . The eccentricity distribution of the thin disk subsample 
suggests that the radial migration played a dominant role in their 
formation and evolution. From the e distribution of the thick disk 
stars it is difficult to conclude which mechanism played the dom- 
inant role (heating, migration, or merger scenario), but the peak 
of the distribution observed at rather lower values than e ~ 0.5 
and the absence of the secondary peak at higher e ~ 0.8 exclude 
the accretion origin. Probably the combination of the two (or 
more) processes such as a heating and migration may describe 
better the observed distribution, just like it e xplains the positi ve 
rotational velocity gradient with metallicity (ICurir et alj|2012h . 

As previously reported by JAdibekvan et al.l (1201 ll) there is 
a family of stars with hig h [g/Fe] values at so lar and super- 
solar metallicities (see also lGazzano et alll2013l) . Although the 
ho-mr stars have intermediate values (even closer to the thin 
disk) of the orbital parameters and space velocities between 
thin and thick disks, they show the same behavior as the thick 
disk stars on the V^ - [Fe/H] and e - [Fe/H] planes. These stars 
share properties of both the thin and thick disk populations. The 
properties of the hamr stars compares well with the me t al-rich 
hig her- [O/H] star s in th e simulation by iLoebman et al.l (1201 ll) 
and iRoskar et al.l (1201 2h . which have intermediate/old age and 
are extreme-migrators (probably < 1-2 kpc). Based on the men- 



tioned results we are inclined to think that the hamr stellar fam- 
ily may have originated from the inner Galactic disk/bulge and 
migrated up to solar neighborhood, although further investiga- 
tions are needed to clarify their exact nature. 
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